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Abstract 

An alternative to the Braginskii decomposition is proposed, one rooted in treating the viscosity 
as a scalar quantity in a coordinate-free representation. With appropriate application to the rate- 
of-shear tensor, one may solve the neoclassical force density equations for its undetermined velocity 
dependence, as well as the radial and poloidal profiles mentioned in [R. W. Johnson, Phys. Plasmas, 
under review], using an improved poloidal expansion. The pseudoplastic behavior of magnetized 
plasma is again obtained, and the high viscosity solution is determined to be physical. A clear 
relationship between confinement mode and viscosity is observed, indicating a physical origin for 
transport barriers, pedestals, and other phenomena. The gyroviscous contribution is found to be 
an effect on the order of one one-thousandth of one percent of the dominant collisional viscosity. 
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I. INTRODUCTION 



Recognizing that quantities which may be written in a coordinate-independent manner 
retain a physical unity that transcends ones choice of representation, but also realizing that 
the magnetization of a plasma introduces potentially complicated effects, we make further 
improvements upon the standard theory of plasma viscosity which was initiated in 

Reference [4] . We isolate the collisional and gyroviscous contributions to the scalar viscosity, 
leaving only its undetermined velocity dependence, the main ion velocities, and the density 
coefficients as the free parameters for which to solve the equations. We also realize the proper 
poloidal expansion for physical quantities in our chosen geometry. We compare the velocity 
dependence of the toroidal and poloidal viscosities for a variety of confinement regimes: L- 
mode, L-mode with Internal Transport Barrier, H-mode, and Quiescent H-mode, and observe 
that the quality of confinement is predicted by the slope of the viscosity-velocity dependence. 
We observe that the internal transport barrier discharge has a toroidal pseudoplasticity 
much like the other L-mode discharge and a poloidal pseudoplasticity more like the H- 
mode shots, suggesting that viscosity is a determining factor of confinement mode. We 
compare the gyroviscous and the collisional viscosity profiles for our collection of discharges 
and find that the collisional viscosity is the dominant effect by a factor on the order of 
O{10^) over the majority of the profile, with slightly more gyroviscous contribution as the 
edge is approached and the various poloidal coefficients approach unity in magnitude. We 
conclude by summarizing and by considering the implications of pseudoplastic plasma on 
fusion tokamak research, as well as on other aspects of plasma physics. 



II. FURTHER DEVELOPMENT 



A. Poloidal Expansion 

We start from the same continuity and force density equations as in Reference jj], except 
for the attention which we will pay to the shear viscosity, and we use the same notation. We 
account for the geometrical effects on physical quantities by first applying the flux surface 
measure factor to the denominator of an expanded quantity, much like the manner in which 
the magnetic field B = B°/(l + ecos6') is expressed, then taking the first order Fourier 
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expansion, eg: 



rij 



(6) = n° (1 + en'j cos 6 + en] sin 9) /{l + e cos 9) , (1) 



where sn^j^ = nj^ has unit normahzation, still neglecting the radial derivatives of the 
poloidal coefficients. With this expansion, the unity, cos, and sin flux surface moments of 
an expanded quantity pick up only their intended term: 

{nj) = —l d9 {1 + e cos 9) Tij (9) =n^. , (2) 
27r Jo 

1 f"^^ 1 _ 

(uj)^ = — (D d9 COS 9 {1 + e cos 9) rij (9) = -en^jrfj , and (3) 

1 1 _ 

{nj)^ = — (p d9 sin 9 {1 + e cos 9) nj{9) = -en%' . (4) 
27r Jo 2 

We note that intermediate results indicate that the choice of expansion has little effect on the 
solution near the core of the plasma but does play a significant role as e = t/Rq increases, 
rendering moot calculations near the edge based on the previous expansion. 



B. Viscosity 

Returning to the definition of the viscous stress tensor (for solenoidal fiows), as given in 
a coordinate-free notation: 

U =7]W , (5) 

where W is the full shear tensor, we note that Braginskii [l| decomposes both the scalar 
viscosity and the shear tensor and then applies certain portions of the decomposed viscosity 
to certain portions of the decomposed shear, to produce the traditional neoclassical viscous 
stress tensor. Without experimental verification, however, this procedure remains to be 
validated. Indeed, unpublished results from the effort behind References 0, 5] indicate 
that the neoclassical model cannot account for the observed velocity profiles in a predictive 
fashion. 

We take the scalar viscosity (V species j) to be the sum of the gyroviscous and the 
collisional contributions, ie 

Vj = Vr + Vf , (6) 

where the gyroviscosity and the collisional viscosity are those of Braginskii (Using 
Shaing's plateau enhanced collision frequency did not produce believable results in this 
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analysis, but more investigation is needed before it can be ruled out entirely.) The shear 
tensor for species j is taken as [sl: 

Wi, ^ d^V^ + d,Vi - ^5.,V ■ , (7) 

where is the Kronecker delta, and is evaluated directly. (Note that the species index can 
be either a sub- or super-script for clarity of notation.) 

Here we depart from Braginskii's model — while Equation ([5]) is correct for a normal fluid, 
as well as for an unmagnetized plasma, the presence of a magnetic field has no small impact 
on the behavior of the plasma. The gyromotion of the ions produces two effects: firstly, all 
the ions on a particular flux-surface execute their gyromotion in lockstep, thus the collisional 
viscosity should be applied only to those terms in the shear tensor without a derivative across 
flux-surfaces; secondly, the ions on adjacent flux-surfaces collide at the gyrofrequency, thus 
the gyroviscosity should be applied only to those terms containing the derivative across 
flux-surfaces. We write our model's stress tensor as 

II, = IIf° + II f = r^f ™Wf ™ + r]f^Wf^ . (8) 

In order to account for the as-yet undetermined velocity dependence, we multiply the vis- 
cosities rij in the a equation by a free parameter, F^, for which we will solve the unity, cos, 
and sin flux surface moment equations. To keep our count of equations and unknowns equal, 
we construct both from a single free parameter via F^ = F° {Pj},. / Pj ■ 



III. APPLICATION AND ANALYSIS 



A. Solution Procedure and Input Profiles 



The astute reader will notice that the toroidal unity moment equation now depends 
explicitly on the poloidal coefficients and will no longer tolerate our previous trick of isolating 
it from the others — we will have to solve the full system of equations at each p. To solve 
our system, we input profiles for density, temperature, poloidal magnetic field (intrinsically 
related to the toroidal current, hence the toroidal electric field as well), toroidal momentum 
injection, and the rotation profiles for at least one species of ion. To analyze shots from DIII- 
D , we retrieve the necessary magnetic geometry and density and temperature profiles from 



EFIT the rotation profiles for Carbon-6 come from GAProfiles 



, and the momentum 



iniection is calculated by NBEAMS l9|. The toroidal electric field is evaluated from the 
plasma current profile using Braginskii's pj resistivity and Wesson's [lyj Coulomb logarithm. 
We neglect the electron and radial flows in the following analysis. Our selection of discharges 
now encompasses a wider variety of confinement regimes, one shot from each of what we 
will nominate as L, L-ITB, H, and QH modes; see Table [T] for the shot numbers and times 



we analyze. All the shots use corrected toroidal velocities 111, and the QH shot also has 

□ 

corrected poloidal velocities [12] ■ The input profiles and radial electric field profiles are given 
in the Appendix. 

B. L Mode Shot 98777 

Let us first analyze shot 98777 in detail. As before, we find two classes of solution, one 
with low viscosity and one with high viscosity, displayed in Figures 1 and 2, but this time 
we will have a criterion for discrimination. As we have applied independent free parameters 
to the toroidal and poloidal components of V • H , the solution is allowed to come to rest 
at independent values for the viscosity associated with flows in the toroidal and poloidal 
directions, but as our theory is rooted in the scalar nature of viscosity, one should demand 
that the physical solution have essentially equal viscosities associated with both of these 
flows. Comparing subfigures (c) and (d) of Figures 1 and 2, we see that the high viscosity 
solution has much better agreement between the toroidal and poloidal viscosities, with the 
velocity dependence accounting for the difference. The viscosities for the deuterium are 
essentially the same line. The low viscosity solution displays a marked difference in the 
viscosities associated with the orthogonal flows, and is discounted as unphysical despite 
some readers' preference for nearly equivalent toroidal velocities between the species. We also 
note the appearance of extreme pseudoplastic or glassy behavior in carbon's high viscosity 
solution, confirming our earlier suspicions That the previous model's high viscosity 
solution is so different than the current one we can only attribute to the previous model's 
improper treatment of viscosity in general (we feel that perhaps a relative - sign should have 
been applied to the toroidal viscosities in which is now provided by the species-dependent 
poloidal coefficients). 

We can investigate the poloidal coefficient profiles for these two solutions. Figures 3 and 
4. Half of the poloidal coefficients have roughly equal magnitude values for the two solu- 
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tions: the density coefficients and the poloidal velocity sin coefficients, with opposite signs 
for the sin coefficients (which are easier to observe in the onhne figures). The variation in 
deuterium's toroidal velocity sin coefficient is accounted for by the variation in the velocity 
profiles. The intriguing difference is in subfigures (d) and (f) — the low viscosity solution has 
nearly equivalent toroidal velocity cos coefficients and poloidal velocity cos coefficients with 
the same shape and a nearly constant factor between the species, while the high viscosity 
solution has deuterium's poloidal velocity cos coefficient of the same shape as its toroidal ve- 
locity cos coefficients and has carbon's poloidal velocity cos coefficient as in the low viscosity 
solution. What these variations signify has yet to be determined. 

We note in passing that the vanishing of the poloidal electric field (when radial electron 
fiows are neglected) established in Reference has consequences beyond the vanishing of 
the poloidal coefficients of the radial electric field — that the electron density sin coefficient 
vanishes is a direct consequence, and empirically we find that the electron density cos coef- 
ficient vanishes as well. We display these vanishing quantities for the high viscosity solution 
for both the new and old models (which stored only one species' determination of Er^^) in 
Figure 5. 

C. Other Confinement Modes 

The velocity and viscosity profiles for the high viscosity solution for the remainder of the 
confinement modes are displayed in Figures 6 through 8. The poloidal coefficient solution 
profiles are given in the Appendix. We note the fairly consistent agreement between the 
toroidal and poloidal viscosities and attribute the differences between the various viscosities 
as a velocity dependent effect driving the values away from a single, common scalar value. 
Pseudoplasticity is observed across the board, with an interesting suppression of glassy 
behavior on the magnetic axis, p = 0. The glassy phenomena is found to occur whenever 
a velocity approaches zero away from the magnetic axis. Glassy peaks are observed in the 
poloidal viscosity for all the shots, and in the toroidal viscosity for shot 99411. These glassy 
peaks near the edge might explain the formation of the edge pedestal and the edge transport 
barrier, and the relatively sudden change in toroidal viscosity in the middle of the profile of 
the L-ITB shot might explain the internal transport barrier. 
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D. Interpretation 

We take the (log of the absolute value of the) toroidal and poloidal viscosities divided 
by the energy density (in other words, Trj/nT = Ft, which is neither the dynamic nor 
the kinematic viscosity, yet appropriate for the current analysis) and plot them against the 
velocities for our collection of confinement modes. This time we combine the deuterium and 
carbon values and separate the toroidal and poloidal values, as displayed in Figure 9. We 
observe both pseudoplastic behavior and the suppression of pseudoplasticity (as evidenced 
by a horizontal relation), particularly in the poloidal viscosity values coming from the inner 
half-radius of p. We take a linear fit of the combined deuterium and carbon data, over 
the whole profile for the toroidal relations and over the outer half-radius for the poloidal 
relations, and label the fitted lines as in the figure. The values for the pseudoplastic slopes 
Sa are given in Table [Tll but as the axes are all on equal scales (relatively, for the velocity 
axes), one may trust ones eye to compare the slopes directly on the figure. We observe that 
the L and L-ITB shots display a greater pseudoplasticity than the H and QH shots, and 
that the values for the toroidal pseudoplasticity seem roughly consistent for both shots of L 
and H mode confinement. The L mode is found to have a greater poloidal pseudoplasticity 
than the others, and we note that the L-ITB has a poloidal pseudoplsaticity more like the 
H and QH modes'. It seems that a measure of the toroidal and poloidal pseudoplasticities 
could be used to define a quality of confinement 

Q = -{S^ + S^) , (9) 

that is on a quantitative scale rather than being a qualitative description, and our values 
for that quantity are in the last column of Table [ITl The lower the value of Q, the better 
the quality of confinement. 

We may also isolate the relative contributions of the gyroviscous and the collisional vis- 
cosities. Figure 10. We observe that there is little variation between the species' toroidal vis- 
cosities and much more variation in the poloidal viscosities, with carbon having the greater 
viscosity. We also observe (note the logarithmic scale) that gyroviscosity is dominated by 
the collisional viscosity by a factor of order (9(10^), hence accounts for one one-thousandth 
of one percent of the total viscosity present. 
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IV. CONCLUSION 



By proposing an alternative to Braginskii's decomposition of the stress tensor, we de- 
termine the pseudoplasticity of magnetized plasma in a fusion tokamak. We observe clear 
relations between pseudoplasticity and confinement mode and suggest a quantitative mea- 
sure of the quality of confinement based thereon. The phenomena of glassy plasma persists 
and is suggested as the physical explanation for a variety of phenomena in magnetized 
plasma, including but not limited to pedestal effects and transport barriers. The imphca- 
tions for other systems of magnetized plasma, both natural and man-made, and for other 
types of plasma, such as the quark-gluon plasma and the gravitational plasma models of 
galaxies, clearly need to be investigated by workers in a variety of fields to determine the 
importance of pseudoplastic viscosity on our understanding of plasma phenomena. We end 
with a conjecture on cosmology — when the early universe went through a phase of being 
electronic plasma, then pseudoplastic viscosity would imply that as flow velocities decreased 
the viscosity would increase until neutralization occurred, impacting the standard comso- 

nnn 

logical model [13|, [l^, 115| in unexpected ways, as well as perhaps seeding the initial density 
fluctuations from which galaxies were born. 
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APPENDIX: ADDITIONAL FIGURES 

1. Input Profiles 

2. Poloidal Coefficients 
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TABLE I: Shot Selection 



Mode Number Time(ms) 

L 98777 1600 

L-ITB 102942 1400 

H 99411 1800 

QH 122338 2750 



Description 
L-mode 

L-mode with Internal Transport Barrier 
H-mode 

Quiescent H-mode 



10 



TABLE II: Pseudoplasticity and Confinement Quality 

Mode Number 5^ Q 

L 98777 -3.73(12) -2.31(11) 6.04 

L-ITB 102942 -3.25(13) -1.93(06) 5.18 

H 99411 -2.28(09) -1.66(04) 3.94 

QH 122338 -2.09(07) -1.32(04) 3.41 
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1. Figure 1. 98777 Velocity and viscosity profiles for the low viscosity solution. 

2. Figure 2. 98777 Velocity and viscosity profiles for the high viscosity solution. 

3. Figure 3. 98777 Poloidal coefficient profiles for the low viscosity solution. 

4. Figure 4. 98777 Poloidal coefficient profiles for the high viscosity solution. 

5. Figure 5. 98777 Vanishing coefficient profiles for both the old and new high viscosity 
solutions. 

6. Figure 6. 99411 Velocity and viscosity profiles for the high viscosity solution. 

7. Figure 7. 102942 Velocity and viscosity profiles for the high viscosity solution. 

8. Figure 8. 122338 Velocity and viscosity profiles for the high viscosity solution. 

9. Figure 9. Viscosity versus velocity for all confinement regimes. 

10. Figure 10. CoUisional viscosity versus gyroviscosity for all confinement regimes. 

11. Figure 11. Radial electric field profiles for all confinement modes. 

12. Figure 12. 98777 Input profiles for poloidal magnetic field, toroidal electric field, 
toroidal momentum injection, C6 toroidal velocity, C6 poloidal velocity, temperature, 
and density. 

13. Figure 13. 99411 Input profiles for poloidal magnetic field, toroidal electric field, 
toroidal momentum injection, C6 toroidal velocity, C6 poloidal velocity, temperature, 
and density. 

14. Figure 14. 102942 Input profiles for poloidal magnetic field, toroidal electric field, 
toroidal momentum injection, C6 toroidal velocity, C6 poloidal velocity, temperature, 
and density. 

15. Figure 15. 122338 Input profiles for poloidal magnetic field, toroidal electric field, 
toroidal momentum injection, C6 toroidal velocity, C6 poloidal velocity, temperature, 
and density. 

16. Figure 16. 99411 Poloidal coefficient profiles for the high viscosity solution. 
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17. Figure 17. 102942 Poloidal coefficient profiles for the high viscosity solution. 

18. Figure 18. 122338 Poloidal coefficient profiles for the high viscosity solution. 



13 





FIG. 1: (Color online). 98777 Velocity and viscosity profiles for the low viscosity solution. 
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FIG. 2: (Color online). 98777 Velocity and viscosity profiles for the high viscosity solution. 
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Shot 098777 @ 1600ms 




FIG. 3: (Color online). 98777 Poloidal coefficient profiles for the low viscosity solution. 



16 



Shot 098777 @ 1600ms 




FIG. 4: (Color online). 98777 Poloidal coefficient profiles for the high viscosity solution. 
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FIG. 5: (Color online). 98777 Vanishing coefficient profiles for both the old and new high viscosity 
solutions. 
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Shot 099411 @ 1800ms 
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FIG. 6: (Color online). 99411 Velocity and viscosity profiles for the high viscosity solution. 
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Shot 102942 (3) 1400ms 
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FIG. 7: (Color online). 102942 Velocity and viscosity profiles for the high viscosity solution. 
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FIG. 8: (Color online). 122338 Velocity and viscosity profiles for the high viscosity solution. 
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FIG. 10: (Color online). Collisional viscosity versus gyroviscosity for all confinement regimes. 
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Shot 098777 @ 1 600ms 
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Shot 099411 @ 1800ms 
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FIG. 11: (Color online). Radial electric field profiles for all confinement modes. 
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Shot 098777(3) 1600ms 




FIG. 12: (Color online). 98777 Input profiles for poloidal magnetic field, toroidal electric field, 
toroidal momentum injection, C6 toroidal velocity, C6 poloidal velocity, temperature, and density. 
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Shot 099411 @ 1800ms 




FIG. 13: (Color online). 99411 Input profiles for poloidal magnetic field, toroidal electric field, 
toroidal momentum injection, C6 toroidal velocity, C6 poloidal velocity, temperature, and density. 
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Shot 1 02942 @ 1400ms 
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FIG. 14: (Color online). 102942 Input profiles for poloidal magnetic field, toroidal electric field, 
toroidal momentum injection, C6 toroidal velocity, C6 poloidal velocity, temperature, and density. 



27 



Shot 122338 @ 2750ms 




FIG. 15: (Color online). 122338 Input profiles for poloidal magnetic field, toroidal electric field, 
toroidal momentum injection, C6 toroidal velocity, C6 poloidal velocity, temperature, and density. 
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Shot 099411 @ 1800ms 




+ D X C 



FIG. 16: (Color online). 99411 Poloidal coefficient profiles for the high viscosity solution. 
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Shot 1 02942 @ 1400ms 
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FIG. 17: (Color online). 102942 Poloidal coefficient profiles for the high viscosity solution. 
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Shot 122338 @ 2750ms 
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FIG. 18: (Color online). 122338 Poloidal coefficient profiles for the high viscosity solution. 
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